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Abstract 

A recent lattice calculation of the spin-dependent structure function g2 is revis- 
ited. It has been recognized that the twist-three operator, which gives rise to c?2, 
mixes non-perturbatively with operators of lower dimensions under renormalization. 
This changes the results substantially. 



*Talk given by G. Schierholz at 'Workshop on Polarized Protons at High Energies - Accelerator 
Challenges and Physics Opportunities^, 17 - 20 May 1999, DESY, Hamburg, Germany. 
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1 Introduction 



The nucleoli's second spin-dependent structure function g 2 is of considerable phenomeno- 
logical interest. The discussion follows the operator product expansion (OPE) In 
leading order of 1/Q 2 , #2 receives contributions from both, twist-two and twist-three op- 
erators. It thus offers the unique possibility of directly assessing higher-twist effects. The 
twist-three operator probes the transverse momentum distribution of the quarks in the 
nucleon, and has no simple parton model interpretation. 

A few years ago we have computed the lowest non-trivial moment of g 2 on the lat- 
tice 0. To convert the lattice numbers to continuum results, the operators entering 
the OPE have to be renormalized. An essential feature of our calculation was that the 
renormalization was done in perturbation theory. It turned out that the twist-three con- 
tribution was the dominant contribution to both, the proton and the neutron structure 
functions. 

Renormalization effects are a major source of systematic error. In this talk we shall 
extend our previous work by employing non-perturbative renormalization. The novel 
feature of this approach is that it introduces mixing with lower- dimensional operators. 

2 OPE and Mixing 

In leading order of 1/Q 2 , and for massless quarks, the moments of g 2 are given by 
2 f dxx n g 2 (x,Q 2 ) = Yl [4% 2 /Q\9(»))di f \») 

(1) 

e{% 2 /Q 2 ,g(»))aH\ri] 



P» n + traces], (2) 

s^Pa)p^ 2 ■•■Pnn + traces] , (3) 

•D^if; -traces, (4) 

■ } ([■■■]) means symmetrization (an- 
tisymmetrization). The operator ([2]) has twist two, whereas the operator (^) has twist 
three. Both, the Wilson coefficients and the operators are renormalized at the scale [i. 
It is assumed that the Wilson coefficients can be computed perturbatively, whereas the 
calculation of the reduced matrix elements a)i and d\i is a problem for the lattice. In 
the following we shall drop the flavor indices, unless they are necessary. 



It -|- i 

f=u,d 



for even n > 2, where [l[] 



and e[^l, S[\ are the Wilson coefficients. Here {■ ■ 
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The lattice calculation splits into two separate tasks. The first task is to compute the 
nucleon matrix elements of the appropriate lattice operators. This was described in detail 
in M. The second task is to renormalize the operators. Generically, 

0(fi) = Z (M 2 ,g(a))0(a). (5) 

As in the continuum, we impose the (MOM-like) renormalization condition 

lTr(g(p)|0( /U )|g(p))[(g(p)|0(a)|g(p))|H _1 = 1, (6) 

where \q(p)) is a quark state of momentum p in Landau gauge. In our earlier work 0, |3| we 
have computed the renormalization constants in perturbation theory to one-loop order. A 
restriction of the perturbative calculation is that it does not allow to incorporate mixing 
with lower-dimensional operators. 

In a recent paper Q we have started a non-perturbative calculation of the renormaliza- 
tion constants associated with the structure functions iq, F 2 and g\. Let us here consider 
the case of the structure function g 2 . We shall restrict ourselves to n = 2. This is the 
lowest moment of g 2 for which the OPE makes a statement. As before [[J, we take o = 2, 
\i x = 1 and /i 2 = 4. We thus need to consider the operators 



<5{ 2 i4} =: {5} (7) 



and 



^f2{l]4} — 202{i4} — Cl{24} — ^4{12} 



-~7 4 D 1 D 2 - ± 74 D 2 Di)j 5 i/; 
=: [5] , (8) 

with O 5 indicating the Euclidean counterpart of O 5 , which belong to the representation 
and t[ 8 \ respectively, of the hypercubic group H(4) ||. The operator @ has dimen- 
sion five and C-parity +. It turns out that there exist two operators of dimension four 
and five, respectively, transforming identically under H(4) and having the same C-parity, 
with which (§) can mix: 

i$(a 13 D 1 -a 43 D 4 S )iP=:0°, (9) 
$( 7l D 3 D 1 - 7l D x D 3 - 74 D 3 D 4 + 14 D 4 =: 0°. (10) 

The operator fllPf ) vanishes in tree approximation between quark states. It therefore 
cannot be included in our framework, and so we discard it. We then remain with 

£> [5] (/i) = Z^(a/j)0 [5] (a) + Z a (an)O a (a). (11) 
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The renormalization constant Z' 5 ' and the mixing coefficient Z a are determined from 
lTr{q(p)\0®(ti)\q(p)) 

lTr(q(p)\0®(v)\q(p)) 
Let us now see how much this effect changes our results for d 2 - 

3 Results and Conclusions 

We work on a 16 3 32 lattice at /3 — 6.0 using quenched Wilson fermions. The lattice 
spacing is a" 1 ~ 1.95 GeV. The calculations are done at three values of the hopping 
parameter k, so that we can extrapolate our results to the chiral limit. We denote the 
reduced matrix elements of the unrenormalized operators and O a by d 2 5 ' and d%, 
respectively. In table |l] we give our results. The numbers for a 2 and d^ (formerly called 
d 2 ) are taken over from 0. These numbers are based on 0(400 — 1000) configurations, 
depending on the value of k. The results for d?, are new, and they are based on O(80) 
configurations only. Hence the relatively large statistical error. 

The calculation of the renormalization constants follows |3J, supplemented by eqs. ([LTD 
- fllcf) . We denote the renormalization constant of the operator by Z^ 5 \ In figs. [[]- 
|3] we plot the renormalization constants as a function of fi 2 . The numbers given refer to 
the MOM scheme, and they have been extrapolated to the chiral limit. We fit and 
Z^ by A + B ln(/i) + C//J,, and Z a by A + C/fi. The result of the fit is shown by the solid 
lines. In order to match our results with the perturbatively known Wilson coefficients 0, 
we need to transform the renormalization constants to the MS scheme. This is done to 
one-loop order H. 



Operator 


0.1515 


0.153 


K 

0.155 


k c = 0.1569 


a? 
af 


0.146(9) 
-0.032(3) 


0.142(9) 
-0.032(3) 


0.138(17) 
-0.044(9) 


0.132(23) 
-0.043(10) 


d f(u) 
df {d) 


-0.097(5) 
0.018(2) 


-0.122(6) 
0.021(2) 


-0.177(14) 
0.032(5) 


-0.206(18) 
0.035(6) 


d a 2 [u) 

4 (d) 


-0.74(12) 
0.136(26) 


-0.89(20) 
0.143(42) 


-1.27(110) 
0.116(119) 


-1.31(69) 
0.137(122) 



Table 1: The unrenormalized, reduced matrix elements 02, d 2 and d?, for u and d quarks 
separately. 



(q(p)\0 [5] (a)\q(p)) 



tree 



— 1 



(12) 



(q(p)\0°(a)\q(p)) 



I tree 



— 1 



(13) 
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fj, 2 [GeV K ] 

Figure 1: The renormalization constant Z^ 5 \ 
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Figure 3: The renormalization constant Z a . 
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Figure 4: The lowest moment of g2(x, Q 2 ) for the proton as a function of Q 2 . 
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We are now ready to state our results for g 2 . At larger values of /i 2 (fi 2 ~ 7GeV 2 ) the 
/i-dependence of and is reasonably well described by the one-loop perturbative 
result (see also 0), so that we may take fi 2 = Q 2 . This choice should keep higher- loop 
corrections to the Wilson coefficients small. In fig. |] we show our new results for the 
proton. The curve marked 'total' gives the total contribution of both, twist-two and 
twist-three matrix elements. The curve marked WW corresponds to the Wandzura- 
Wilczek approximation J7| , in which the contribution from twist-three matrix elements is 
discarded. The dashed lines indicate the errors on the curves. The errors come almost 
completely from the error of the nucleon matrix elements. 

Unfortunately, our errors on d 2 are yet too large to reach a definite conclusion. But 
we may say already that the effect of mixing is large, and of the same magnitude as the 
twist-three contribution itself. At Q 2 = 5 GeV 2 we now find for the proton 

/ dxx 2 g 2 (x,Q 2 ) = 0.007(15), (14) 
Jo 

and for the neutron we obtain 

f dxx 2 g 2 (x,Q 2 ) = 0.006(7), (15) 
Jo 

which is not inconsistent with the experimental values || . In case of the proton it appears 
that d 2 vanishes at larger values of Q 2 , leaving us with the Wandzura-Wilczek result. 

We hope to return with more precise numbers in the near future. 
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